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SUMMARY 



Elenientary basic solutions of the equations of 
motion of a compressible fluid In the hodograph variables 
are developed and used to provide a basis for comparison, 
in the form of velocity correction formulas, of corre- 
sponding compressible and incompressible flows. ^The 
known approximate results of Chaplygin, von Karman and 
Tsien, Temple and Yarwood, and Prandtl and Glauert are 
unified by m_eans of the analysis of the present paper. 
Two new types of approximations, obtained from the basic 
solutions, are introduced; they possess certain desirable 
features of the other approximations and appear prefera- 
ble as a basis for extrapolation into the range of high 
stream Mach numbers and large disturbances to the main 
stream. Tables and figures giving velocity and pressure- , 
coefficient correction factors are included in order to 
facilitate the practical application of the results. 



INTRODUCTION 



The present paper is concerned with a theoretical 
study of the hydro dynamical equations of a perfect com- 
pressible fluid in two dimensions, in which the so-called 
hodograph variables are used as the independent variables. 
It is hoped to achieve herein a unification of the 
present-day results obtained in this field and also to 
provide a working basis for further developments. The 
earliest contributors to the hodograph m.ethod for treating 
compressible fluids were Molenbroek (reference 1) and 
Chaplygin (reference 2) . The remarkable work of 
Chaplygin on gas jets ajjpeared in Russian in 1904 but re- 
mained relatively unnoticed. In recent years contribu- 
tions to the hodograph method have been made chiefly by 
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Demtchenko (reference 3), von Karman (reference 4), Tslen 
(reference 5), P^lngleb (reference 6), and Temp?Le and 
Yarwood (reference 7) . ■ 

The chief reason,, and perhaps the only reason, for 



preferring th 
coordinates i 
hodograph var 
achieved, how 
conditions an 
siraplif icatic 
the less makes 
as long, as it 



e hodograph variables to the physical plane 
s that the equations of motion in the 
lables are linear. This simplification is 
ever, at the cost of more difficult bo -on dary 
d at a loss of physical insight. The great 
n in the mathematics due to linearity never- 
it desirable to pursue this line of attack 
aT.>pears profitable to do so. 



The mathematics for handling the flov^ equations re- 
ceived a substantial impetus by the work of Bers axid 
Gelbart (reference 8), who developed a new function theory 
analogous to ordinary analytic function theory. The 
present paper utilizes the methods of this new function 
theory to develop certain functions essential to the 



c omp r e s s ib 1 e - f 1 o w pr ob 1 em^ . 



.s of historical Interest 



that Ideas similar to those of Bers and Ge3_bart were ex- 
plored by the renovmed mathematician HUbert (reference 9) 
In the early part of tliis century but do not appear to 
have been further developed at the time. 



The m.aterial 
into two parts. 



to be treated is conveniently separated 
In part I, the present paper, basic 
particular solutions of the hodograph flov/ equations are 
developed and employed in unifying and extending the re- 
sults obtained by Chaplygin, von K^rman, and Temple and 
Yarwood, The results obtained in part I are of immediate 
practical application and are given in the form of tables 
and graphs' of velocity and pressure-coefficient correction 
factors. In part IT, which will appear later, general 
particular solutions of the hodograph flow equations are 
developed and discussed. The material in part II, it 
Is hoped, will lead to a method for handling the actual 
bomdary problem of the flow of a com.pressible fluid past 
a prescribed body. 



ANALYSIS 

Plow Equations of an Incompressible Fluid 

It is v;ell known that the relations between the 
velocity potential 0 and the stream^ function >j/ for 
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the steady irrotational two-dimensional motion of a 
perfect incompressilDle fluid are 

5x dy ' 



5y " 6x 



(1) 



These equations are the Caixchy-Riemann equations and 
therefore ^ + i\J; is an analytic function f(z) of the 
complex variable z - x + iy. 

• The complex velocity or reflected velocity vector 
u - iv is obtained from the complex potential f(z) by 
differentiation. Thus, 

df (z) 



qe-i9 

^-i(9+llogq) 



(2) 



where q is the magnitude of the velocity vector and 8 
is the angle the vector makes with, the positive direction 
of the X-axis. 

The variables 9 and q are sometimes referred to. 
as- the hodograph variables.^* The flow equations in the 
variables 9 and q can be readily derived by intro- 
ducing 9 + i log q as the independent complex variable 
in -place of x + iy. Then, in analogy with equation (1), 



d log q 



60 



6 log q 69 



(3) 



mmmsmmm 
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or 



de 



= q 



6i 



dq 

1 d± 

q 69 



(4) 



These equations are knoTO as the hodograph equations for 
the flov/ of an incompressible fluid. 



Flow Equations of a Compressible ir'luid 

The equations corresponding to equation (1) are, for 
a compressible fluid, 



dx p 6y 
6^ ^ ^fo 6\j/ 



(5) 



where p Is the density of the fluid at any point (x,y) 
and is a constant density, which for convenience is 

referred to a stagnation j)olnt» 

A short way to derive the hodograph equations for a 
compressible fluid, attributed to Molenbroek, is as fol- 
lov/s : 

According to equation (5). with u ^ and v=t^ 

ox oy 

d)5 + i ~ dvj/ = (u dx + v dy) + i(-v dx + u dy) 

(u - iv) (dx + i. dy) 
= qe-^Q dz 



or 



dz = ^ e^^(di?^ + 1 ^ (6) 



• • • 

• •• 
• 



• • • 
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It follows from equation (6), by considering 6 and q 
as Independent variables, that 



and 



dz 
69 



6z 



q ® \de ^ P 68/ 



1 eie 
q 



\dq " P 6q/ 



Then, by assuming that p is a function of only q 
(equivalent to assuming that the pressure is a function 
of only the density), 



6q 69 



and 



69 6q 



= e 



19 



J:_ 6^ . d(Po/Pq) d£\ 1 i9/V£. , . Po 
_"q2 6e ^ ^q 6ej ^ ® \6q 69 P 6q 69/ 



q \6q P 6qy q. \6Q 6' 



q ^ P 69 dq 



Since, by continuity, these two expressions are identical, 
it follows that 



q ^ V6q P 



q2 69 dq 69_ 



Hence, by equating real and imaginary parts. 



69 P 6q 

M = ^(Pq/p^) 6j/ 
6q ^ dq 69 



1 



> 



(7) 
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These are the hodograph equations, first obtained by 
Molenbroek, for the fiov/ of a compressible fluid and are 
independent of the form of the pressure-density relation. 
It is observed that, when P = Pq - Constant, equa- 
tions (7) reduce to equations (4), Equations (7)^ in 
contrast v^ith equations (5), are linear in the dependent 
variables . 

Bernoulli's Equation and Equation of State ^ 

In the present section there is listed a collection 
of formulas and definitions necessary in the analysis. 

Bernoulli's equation for a compressible fluid is 

o 

where 

p static pressure in fluid 

p_ static pressure at stagnation point (q =^ 0) 

o 

p density of fluid 

q magnitude of velocity of fluid 

The adiabatic relation between the pressure and the 
density is 




where 

Y adiabatic index (approx* 1,4 for air) 

Pq . density of fluid at stagnation point (q = 0) 

The local velocity of sound a is obtained from 

dp 
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= Y ^ 



(10) 



Prom Bernoulli's equation (8) and from equations (9) 
and (10), the follov'ing relations maybe obtained; 



.2 _ 



ao^ - -i (y - i)q' 



P = P 



oj .2 

L. 



P, 



1 - (Y - 1) 



q3 


Y-1 






q2 


Y 


Y~l 


R 2 





(11) 



^J^r■^,eve a^^ is the velocity of sound at stagnation point 
(q = 0) 

From equations (11), for Y > 1^ a iriaximum velocity 
q - q^^^ is obtained for the limiting conditions 



p = p = a = 0. 



L hus , 



where 



a. 



'^m V - 1 0 
^ 2pa 2 



(12) 



" Y - 1 

The fundamental nondirnensional speed variable, in 
general, is q/a^ but It is found useful in the analysis 
to employ a nondirnensional speed variable t defined as 



•m 



(13) 



For Y > 1^ the range of the variable t is 0 ^ t ^ 1. 
The value t = 0 has a dual meaning; t o in the case 



• • • • < 

• • • • • • • 



• • • 



• • ••• 



• • • • 
• • • 



I • • 
I • • 
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of a coinpresslble fluid corresponds to a stafrnation point 
(q = 0), or T = 0 may mean the limiting case of an in- 
compressible fluid {Rq = . 

With the definitions of t' and p, equations (11) 
become 



1/2 



a = Bod - t) 
p = Pq(1 - t)P I> 



(14) 



The local Bilach number M = q/a may be expressed in 
terms of the speed variable' t in the follov/ing way: 



m2 = ^ ^ ^ 

a^2 a 2 
o 



1 - T 



(15) 



or, by solving for t in terms of M, 



T = 



2p + 



(16) 



The value of t for v;hich the local velocity of 
the fluid equals the local velocity of sound (M = 1) is 
given by 



2p + 1 



(17) 



In the case of uniform flow past a fixed boundary, 
the pressure coefficient is defined as 



'P,Mi 



= ^ " Pi 

2 Pl^l 
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viThere the siibscrlpt 1 refers to the vmdlsturbed stream, 
The pressure coefficient for the incompressible case 
in = 0) is 



(18a) 



The pressure coefficient for the compressible case is 





For q = qg (sonic) , 



<'-l + 



Y 



+ (y - i)M]_^!y-i 



Y + 1 



J 



> 

J 



(ISc) 



For ^ - (vacuum). 



(18d) 



Basic Solutions of Hodograph Equations 

Consider the incompressible case represented by 
equations (3) or (4). It is clear that 0 = Q and 
\j/ = log q satisfy these equations * In fact, any con- 
vergent power series in w = 9 + i log q represents an 
analytic function of vi/hich the real and imaginary parts 
satisfy equations (3) or (4). The class of analytic 
functions in w (and the concept of analytic continua- 
tion) then yields all the particular solutions of these 
equations « 



The particular solution = 9 +1 log q can be 
obtained by means of an integration that is instructive 
in the generalisation to the compressible case. It is 
well known that 



• • • • 

• • • • • • 



• • • 



• • • • • 



• • • • 
• • • 
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P(w) =jf(w) dw 
can "be represented as the sum of tv/o line integrals 

P(w) = ^(P d8 - Q dlogq) + xj (Q de + P d log q) 



where 



■(w) = P + IQ 



Thus, given a pair of functions P and Q that. satisfy 
equations (3) or (4), thin process yields another pair 
of solutions, namely, the real and the imaginary parts 
of P(w). For eixample, if P = 1 and Q"= 0, 



P(w) =: w 6 + i log q 
Again, if P = 0 and Q, - 1, 

P(w) = iv/ = -log q + ie 



(19) 



(20) 



The physical interpretation of equations (19) and (20), 
considered as flow patterns, is of some interest in con.- 
nectipn with later developrients , It is clear that 
equations (19) and (20) represent a vortex and a source 
located at the origin, respectiveljr. 

The generalization to the compressible case of the 
foregoing elementary results was accomplished by Bars 
and Gelbart (reference 8) by means of simple yet fertile 
ideas. Bers and Gelbart treat equations, of the form 



b0 
•58 



6q 



-X2(q) II 



(21) 



and show as is readily verified that^ if P and Q are 
a pair of solutions ^ the real and imaginary parts of the 
following 3un of line integrals 



• • • • ••• • 

••• • 

••••• • ••• ••••••••• 

• • • • • • 
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y^[p cl9 - \^{q) Q dq] + Q dO + —i-^ P dqj (22) 

are also solutions of equations (21) . 

In particular, corresponding to the pair of solu- 
tions P = 1 and Q = 0, there is obtained 



and, for P = 0 and Q = 1, 



.1 



IW = i|e + ij\2(q) dqj 



(24) 



By repeated application of the process of integration, 
indicated by expression (22), a general set of particular 
solutions of equations (21) may be obtained. These 
particular solutions are discussed in part II; in the 
present paper, only the solutions given by equations (25) 
and (24) are needed. 

The general hodograph equations (7) are of the form 
of equations (21) with 

and 

d(Po/pq) 



Xniq) = -q 



dq 



For the rest of this paper, the adiabatic pressure-density 
relation (9) is used. 3y means of equations (9) and (14) 
and the relation 



dq q 
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obtained from the differential form of Bernoulli's equa- 
tion (8), it follows that 



(1 - T)P 

_ 1 - (2(3 + 1)T 
q(l - T)P+l" 



(25) 



The evaluation of the integrals in equations (23) and (24) 
is made unique by requiring that the results reduce to the 
incompressible case when the speed of sound la infinite 
(that is, when t = 0) . Then, 

L = J(l - T)P^ 

= log q + f (t) (26) 



where 



o 



and 



J ^ ^)(3+i q 



= log q + g(T) 



(27) 



where 



g(T) =1 

L (l-T)H 



1 - (2 p + 1) t 
+1 



- 1 



en 
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and it is observed that the functions f(T) and g(T) 
vanish for t = 0, 

Equations (25) and (24) can be written in the form 

?/ = 9 + IL 

and 

iW = i(e + iL) 

It is Important to note that, in the incompressible case^ 

W and iW reduce to and iw, since L and L re- 

duce to log q. Thus, there are in the compressible 
case two basic functions L and L corresponding to the 
one function log q in the incompressible case* It is 
of interest to mention that the functions V/ and VN , 
considered as flov/ patterns in a compressible fluid, can 
again be interpreted as a vortex and a source - 



Evaluation of Functions f(T) and gd) 

for Various Values of p 

In general, the integrals in equations (25) and (27) 
representing the functions f(t) and g(T) are expres- 
sible by infinite series. For the important case of air, 
hov/ever, with the adiabatic index r put equal to 1.4 
instead of the usual value 1^408, these functions can be 
obtained in closed forms. Thus, v/ith p = 2.5, 



f(T) = 




(1 
+ (1 



- T)5/2 ^ i (1 - T)3/2 



- 7)1/2-. 23 1 + (1 - T)^/2 

} 15 2 



(28) 



• • • • • ••• 

«• • ••• ••• ••• 

••••••••• 

••• • ••• • • • • • 

• « • • • ••• •• 
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and 



, , 1 /"T 1 -, 6T 



T 



1 ^ 1 L 

- T)^ - (1 - T)"^ 



(1 - T)"/'^ - (1 - T)" 

1_ 1 - log ^)'^' (29) 

(1 - T)l/2 3 ^ 2 



Table 1 contains values of f(T) and g(T)^ and fig-- 
ure 1(a) shows, these functions plotted against r. ' Ob- 
serve that f(T) and g(.T) are well-behaved fiinctions 

in the range 0 ^ r < 1. in figure 1(b), these functions 
are plotted against the local Mach number M in the 
practical speed range. 

Other interesting cases for which the functions f (t) 
and g(T) can be expressed in closed forms are y - ^> 
Y - 2^ Y = 3/2, and y - -1. For r - ^ (p = 0, 
a = incompressible case), 

f(T) = g(T) 0 



For Y = 2 (p 1), 



^(t) ^ 



g(T) = 1 - Y'^^ " i (1 - t) 



For Y 



1 (p 



= 2), 



f(T) = 



X ^ l t2 

T + 4 T 



> \ 1 1 1 + T 1 ^ . 

=2-2 (1 > T)2 " 2 (1 



- T) 
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For Y = -1 (p = 4)' 



f(T) = g(T) = - log 



1 + (1 - T) 



For the Isothermal case Y =^ 1 (P the velocity 

of sound a = ao = Constant and the functions f and g 

are obtained infl-nite series in the ratio q/aQ. Thus, 
in the limit p-r->op^ ■ ' ' 



■ Lim 1 



1 ^ 



28a, 



- 1 



2 Jd(q/ao) 

qAo 



- 1 



<^(qAo) 
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For arbitrary values of y (or .p) the expressions 
for f(T) and g(T), obtained v^ith the aid of the bi- 
nomial expansion, are 



n=l 



4 a 2 32 ^^W2 
^o vo 



and 




The significant feature of this general result is that, 
if powers of q/ag higher than the third are neglected, 

f(T) = g(T) -i • (30) 

^O 



and does not involve explicitly the adiabatic index y. 
This circumstance underlies the present-day approximate 
methods for obtaining velocity and pressure-coefficient 
correction factors; in the following sections, this point 
is brought out more clearly. 
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Application of Basic B'lanctions L and L 

In this section, the basic functions L and L are 
employed to set up relations .between velocities in 
^^corresponding^* compressible and incompressible flows. 
These relations are of the nature of "stretching factors" 
or velocity correction^ formulas and contain the results 
of Chaplygin, von Karman, Temple and Yarwood, and Glauert 
and Prandtl. It is important to recognize at the outset 
that no single velocity correction formula can represent 
in an exact v/ay the correspondence of flow patterns past 
a prescribed body in a compressible and an incompressible 
fluid, A single velocity correction formula is actually 
feasible in only two cases: (1) The stream Mach number 
is small (even though the disturbance to the main stream 
due to the presence of the body may be large) so that the 
compressible-flow pattern differs only slightly from the 
incompressible-flow rjattern or (2) the disturbance to the 
main stream is vani shingly small (even though the stream 
Mach number may be high) so that the effect^'of the shape 
of the solid boundary is small. The various velocity 
correction formulas discussed in the present paper differ 
essentially only in the degree to which the requirements 
of these two cases are satisfied. Despite their limi-. 
tations, single velocity correction formulas are extra- 
polated, in viev; of the lack of more rigorous solutions, 
into the range of large disturbances to the main stream 
and high Mach numbers. This extrapolation can be 
justified by further theoretical investigations and by 
comparison v/ith experimental results. 

Consider again the corresponding pairs of functions 

w = 9 + i log q 
Vv' = e + iL 

and 

iw = 1(9 + i log q) 
if = 1(9 + iL) 



It has previously been noted that the pairs of functions 
in equations (31) and (52) denote respectively a vortex 
and a source in an incompressible and a compressible 



J 



(31) 



(32) 
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fluid. Each pair of functions can be employed to def5,ne 
a correspondence of flow patterns in v^hich corresponding 
points are identified by the same values Thus, 
in the case of the vortex (equations (51)), 



01 = = 9 
^i ^c ~ ^1 = L 



where the subscripts i and c refer to the incompres- 
sible and to the compressible case, respectively. It 
follows that 



fir) 



(33) 



Similarly, in the case of the source (equations (32)), 



0^=0^= -log = -L 

and 



qi = 



= .^0 



e 



g(T) 



(54) 



At the end of the preceding section it was pointed out 
that, to a first approximation^ the functions f(T) and 
g(T) are equal. This fact implies that, to a first 
approximation, a single velocity correction formula is 
feasible. The assiomption is nov; made that either equa- 
tion (33) or equation (34) can be adopted to provide a 
correspondence of flow patterns in the case of uniform 
flow past a body in an incompressible and a compressible 
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fluid. Vv'ith the undisturbed streams as convenient refer- 
ences, the following nondlmenslonal forms of equations (33) 
and (34) can be written: 



and 



O ^ 1 ' 



(36) 



e 



where the subscript 1 refers to the undisturbed stream. 

The use of the undisturbed stream as reference in the 
nondimensional form of the velocity correction forraula 
was introduced by Tsien in reference 5, v/here also the 
details of the von Karman approximation are developedo 
It is shomTi in the follov/ing section that either of 
equations (35) or (36) contains the results of Chaplygin, 
von Karma!n, and Temple and Yarwood. As has been previ- 
ously pointed out, the concept of a single velocity cor- 
rection formula is feasible in only two cases, namely, 
small stream Mach numbers and vanishingly small disturb- 
ances to the main stream. It is desirable then to seek 
a single velocit^r correction formula that com.bines the 
features of these two cases. Prom this point of view, 
equation (35) or equation (36) is not the best choice • 
A better choice of a single velocity correction formula 
appears to be the following combination of equations (35) 
and (36), based on the arithmetic mean of f(T) and g(T): 

|Cf(T)+g(T)J 

/JL\ . ^ (37) 

V^lA Ufir,)^gir,)\ 



e2 



In a later section, still another combination refer- 
red to as "the geometric-m.ean type of approxim.ation" is 
introduced; in the section dealing with the Glauert-Prandtl 
approximation, certain features of the foregoing arithmetic- 
mean type of approximation and of the geometric-mean type 
are discussed. 
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At this point it is desirable to discuss the practical 
application of equation (37). According to equation (16), 

2p + 



and 




(38) 



Equation (37) then yields, for a given set of values of 
the stream Mach number M]_ and the local Mach number M, 
a value for the ratio {^^^l)± local velocity q 

and the stream velocity q^^ in an incompressible fluid. 
Table 2 shows corresponding values of (jl/^l)^ 
(^^/O-l^l for various values of the stream Mach niomber Mj 

with Y ='1.4 (P = 2.5). This tabulation is performed, 
for the purpose of comparison, for the three cases repre- 
sented by equations (35), (36), and (37). Values of 

(q/qi)3_, (^/^^0.)c^ ^"'^ (q/qi) >^ obtained from equa- 
tions • (37) and (38), are plotted against the local Mach 
number M in figure 2 for various values of the stream 
Mach number Mi* Table 2 als'o shows values of the 

pressure coefficients Cp^o ^Vp^l calculated by 

equations (18a) and (18b) for these corresponding values 
of ^^^^ C^^^^Oc' Figure 3 shcv/s "the curves 

of pressure coefficients corresponding to the curves of 
velocities of figure 2. Useful cross plots of the 
curves in figure 3 are shown in figure 4, in v/hich Cp^iji^ 

is plotted against Mi for various values of Cp^^. In 
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addition, curves are shown ' In .figure 4 for (Cp^^i^g 
and (^p,Mi^j^ calculated by equations (18c) and (18d) , 
respectively. The curve for (Cp^Mi)g corresponds to 
the sonic value M = 1 or t = Tg = l/6 and in effect 
divides the region of flow into a subsonic and a super- . 
sonic part. The curve of (Cr. ) corresponds to the 

maximum value M = or T = 1 and represents the outer 
limit of the supersonic, region (or a perfect vacuum) • 
In order to exhibit the main differences between the 
various correction formulas (35), (36), and (37), the 
ratios of the sonic values (Cp^Mi)^ the correspond- 

ing incompressible values Cp^o are plotted against the 
stream Mach nuraber in flgiire 5o 

Observe in figure 2 that the fq/q^ ), -curves have 

maximum points. This fact means that the value of 
(q/^l^^ associated with a value of i^^/^lji 

unique. Analytically, the criterion for the maximum 
point is equivalent to 




(39) 



or, from velocity correction formula (37), 
(1 ^ T)^P'*^ - (2p + 1)T + 1 = 0 



For p = 2.5 this equation has only one positive root, 
T « 5/24 or M 1.15. It is interesting to note that 
velocity correction form.ula (36) yields as 'the criterion 
for the maximum point 

1 - (2p + 1)t = 0 



The root of this equation is t = 73 = ■ , and, for 

^ p + 1 ■ ■ ' 

P = 2.5, is T = 1/6 or M = la Velocity correction 
formula (35) yields no maximum value of t or M. 
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Meaning can be given to the value t = l/6 (M = 1) 
in the case of equation (Jm-) '^^ith reference to the origi- 
nal Interpretation of the flo\v pattern as that of a 
source. It can be shovm that the acceleration 

(q along a streamline is infinite at all points 

for which the local Mach number is unity (r = l/6) and 
that a flow discontinuity exists there. In the case 
of the vortex flov/ pattern (equation (35) no flow 
discontinuity occurs for M< <» . The velocity correc- 
tion formula (57) suggests a ^'limiting" value M ^ I.I5 
for a spiral flov;, since equation (59) is analogous 
to a condition of infinite acceleration, Ttms, the 
existence of a mixed subsonic and supersonic region of 
flow ^-vithout diacontinuities is indicated. Since the 
occurrence of this liiiiiting valu.e of M is a consequence 
of the simple form assumed for the velocity correction 
formula, no undue significance should be attached to any 
particular value at the present time. 



The Chaplygin Approximation 

From the point of view of the present paper, 
Chaplygin 's approximation for subsonic speeds assumes a 
simple and lucid form. Chaplygin introduces in place 
of q a new independent speed variable equivalent 
to the quantity given on the right-hand side of equa- 
tion (55) J nam.ely, 

T) = qe^^''^^ 

The hodograph flow equations (7) then assume the form 

1 

t 

> 



68 ' 6t) 



(Ii-0) 



6ti 



= -.P(T) 



61 



J 



v\rhere 
P(t ) 



1 - (2P + 1)t 
(1 - tVp-^I 
1 - p{2p + 1)t2 



I p(2p + 1) (2p + 2)t5 - ... 



• •• 

• •• 

• • 

• • • 

• •• 



• • • 

• • • « 
•••••• 



• • « 

• • • i 
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Values of the function ?(t), for several values of Y 
(or (3), are given In table 3 and are plotted in figure 6 
attains t the local Mach number Mo Chaplygin noted that, 
in the case of air (p = 2.5), F(t) differs but little 
from unity over about one-half the subsonic range 
0 = T = l/6. His approximation in the range of lov; sub- 
sonic speeds consists in neglecting powers of t higher 
than the first or in replacing F(t) by unity. Equa- 
tions (40) can then be written in the Cauchy-Riemann 
form 

M = 

69 d log T) 
6 log T} ^9 

and + i\|/ therefore is an analytic function of the com- 
plex variable 9 + i log y]. Chaplygin 's approximation 
thus leads to the velocity correction formula 



1 - I T 



v«here powers of t higher than the first are neglected 
throup-hout. The use of equation (34) instead of equa- 
tion T33) also leads to this result to the same order of 
approximation. 



The von Karman Approximation 

Von Karman 's approximation corresponds to the case 
Y = -1 (or p = . It follows at once from the 
integral expressions for f (t) and g(T) given by equa- 
tions (26) and (27), respectively, that for this case 



f(T) = g(T) = - log 



1 + (1 ~ T) 



1/2 
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or, with the use of equation (16), 



f(T) = g(T) = - log 



1 + 



(i - M2y/2 



This function, plotted against M, is included in. fig- 
ure 1(b). Corresponding to equations (35) and (36), 
there is a single equation 



51 /i V^l/c I +(l - T )l/2 



Ml 



Replacing t by Tt/-^\ and t-i by — according 

to equation (16) yields 



1 + (^1 - M-L^^' 



1/2 



(l - . 



1 - Mi^ + Mi^^^^ 



2 



a/2 



(42) 



Then, by solving for {^/'ix)^ terms of (^/'^D 
the stream Mach number M]_, 



vfhere 
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The pressure coefficient ^p,Mi' expressed in terms 
of the incompressible pressure coefficient ^'p,o> 
easily obtained from the general formula (18b) by putting 
Y = -^l and making use of equations (43) and (18a). 
Thus , 



(l - Mi2) 



1/2 



(44) 



1 + (l 



- Mi2) 



.1/2 2 



Observe that for this case the function P(t) intro- 
duced b;/ Chaplygin and given in equation (40) is exactly 
equal to unity. ^From the point of view of the present 
paper then, von Karman^s approximation appears to be 
equivalent to that of Chaplygin, who approximates F(t) 
by unity. It follows that the range of validity of 
von Kantian ^-^i approximation and that of Chaplygin, in a 
strict sense, coincide. Furthermore, it is pointed out 
that the von Karmf^n approximation does not permit a 
supersonic region. Von Karm^an's choice of y - -1 i^^s 
the advantage, hov/ever, of yielding simple explicit ex- 
pressions for (q/q^^^ in terms of ^q/q^^^ and for 

Cp^j,/T^ in terms of G^^q. Several values of Cp^j/[^ 

calculated by equation (44) are included in figure 4. 

For the purpose of comparison with the other approximations, 

there is plotted in figure 5 the ratio of (^P,Mi)g 

Cp Q against the stream Mach nuraber M]_ in the case of 
von Karm.an's approximation. The values of C^ ^ are 

obtained with the use of velocity correction formula (42) 
for the local Mach number M = 1, but the values of 
C^P ^*^l)o calculated \vith y = lo4 6 



The Temple-Yarwood Approximation 

The functions 0 and if related by the first-order 
simultaneous equations (21) separately satisfy the second- 
order equations 



• • • • • • • ••• ••• 

• •••••••• 
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_J.-.M!= 0 



6^^!/ ^ 1 



~ i I q / 



= 0 



In teriiis of the nondimensional speed variable t and with 
the values of X;j_(q) and \2(<i) adiabatic case 

given by equations (25), these equations take the form 



\ 



1 (1 " t)P rf^ ^ b r_rji -_t)P_"^3. 5^^" 



= 0 



S (^6) 



1 1 - (2p. + 1)t , .6' 

k /, A'>-2 6t 

t(1 - T ) Cc 



(1 - t) 



p 01 



= 0 



Pormal solutions of these equations wei'e given by Chaplygin 
in the form of tv/o infinite series 



>i! = B8 +2_ Bja sindae + e^) 



> (ii-7) 



where the functions N'^'mtT) and 0y,^i'^ ) are obtained from 
hypergeorietric series and B, Evj, " and €^ are arbitrary 



constants , 



A disadvantage of the formal solution, as rernarked 
by Temple and Yarwood, is that it is unsuitable for 



• • • 

• •• 



NACA ACR No. L4C24 



CONFIDENTIAL 



27 



numerical computation "because the hypergeornetric functions 
Involved are complicated and are not tabulated. Temple 
and Yarv^^ood therefore looked for approximations that are 
of practical value In calculations of compressible flows. 
By means of a skillful analysis, they found such approxi- 



mat ions and showed that the 
are of the type 



simplest forms for \]/^ and 



m 



nm 



4(T) s [S(T)J 



4(t) log g(T) 



(48) 



where r\{y) and , Independent of the index m, are 



(49) 



Significantly, from the point of view of the analysis 

of the present paper, the functions rj and approxi- 

functions 



m.ated by ^1 - "^y 



are none other than the 



defined on the right-hand sides of equations (35) and 
(34) • The approximation of Temple and Yarv/ood then 
leads to the same velocity correction relation as was 
obtained by m.eans of Chaplygin's approximation (equa- 
tion (41)) . 

The velocity and pressure-coefficient correction 
formulas obtained by Temple and Yarwood are more involved 
than the explicit expressions (45) and (44) obtained by 



von Kc 



rmAn. 
2 



Replacing 
thus yields 



T in equation (41) by 



:/iL\ 



1 - ^ 



1 - 



5 
4 



(50) 
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v/here 



^1 = 



5 + M' 



The solution of this cubic equation for (<l/q-|^^ is 



(tl - {AX <^ - i -) 



COS ^ (rr + a) 
""^ covS a 



(51) 



where 



CCS a = 



and 0 < a ^ ~ 



The pressure coefficient C> 



is 



then calculated by equation (18b). Some values of the 



pressure coefficient C 



calculated with the aid of 



equation (51) are shown in figiire 4; a curve of 
i^9,Ki^^^^p,o plotted against is included in fig- 

ure 5. It is remarked that^ with the use of equa- 
tion (39), the velocity correction formula (50) yields a 
limiting value M =: 1.35. 



Approxiination Based on C-eometric Mean of dL and dL 

Vv'ithout going into its deep significance in the 
present paper, it Isof interest to introduce another 
function related to L and L and to the general 
particular solutions. This function, which like L 
and L reduces to log q for t = 0, is defined by 



H(T) = J (dL df )^/^ 



(52) 



It is remarked that H(t) is closely related to a 
function K(t) employed by Temple and Yarwood (refer- 
ence 7) in the determination of their approximation. 
In the next section, it xiill be seen that the function 



• • • • 
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H(t) plays an important role in connection with the 
P r andt 1 - G 1 a ue r t app r o x ima t i on t 

Prom equations (26) and (27), 



dL = 



^d,= (l-r)3f 



and 



(1 - T)P-^1 



Then, 



(dL dL) 



\l/2 



i 1 - 



p + 1)t1^/^ 



1 - T .J 



dq 

q 



(53) 



and, from equation (52), 



H(t) = log q + h(T) 



(54) 



where 



h(T) 



1 
2 




1 1 ai 
J ' 



The function h(T) can be obtained in a closed form for 
any value of y (or p) and is 



h(T) = -log 



: ■ ^ , (55a) 



2(1 - V-isj 



1^ 
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where Tg = l ^^'^ v/here this expression is valid 

in the i^ubsonic range 0 ^ t ^ Tg, V.'ith t replaced 
by — and 0 ^ M ^ 1, the expression for h(T) 

becomes 

hfT) - lo. 1 ^ (1 - m2)'1/2 ^ _ 1.1,2)1/^ 

h^.T) - -log --7=r- log == 

+ — ^== — log :; j= (55b) 

2Vt^ 1 + 



It is observed that, for the supersonic region 

Tg ^ T ^ 1 or M > 1, H(t) as defined by equation (52) 

becomes a complex functionj but, for present purposes, 
only the real function of the subsonic range is utilized. 

The function H(t) may be utilized to obtain a 
velocity correction formula in the same manner as the 
functions L(t) and L(t), Thus, analogous to equa- 
tion (35), (36), or (37), 




(56) 



It is instructive to compare equation (56) with the 
approximation given by equation (37). Equation (37) may 
be written as 



f-J(dL+dL) 
/_q_\ _ e^-^ 



(dL+dL) 



T=:T' 
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and equation (56) may be written as 

/(dL dC)^/^ 



e 



Thus, the pov/er of the exponential is in one case the 

*^ 

integral of the arithmetic mean ^^i-i-JlIi and in the 

other case the integral of the geoirietric mean 

(dL dL)-^/^. Table 1 shows values of the furtctlons 

lllLj^.3lll and h(T) in the case of air (y = 1-4, 

p = 2.5, and T3 = 1/6) and figures 1(a) and 1(b) show 

these functions plotted against t and M, respec- 
tively. Observe that these functions, and consequ.ently 
the velocity correction formulas (S7) and (56), differ 
only slightly in the subsonic range 0 < M < 1. Pig- " 
ure 5 exhibits graphically a comparison of the velocity 
correction formulas (37) and (56) for M = 1. The 
limiting value of M (defined by equation (39)) is 
M =1 in the case of equation (56) as compared with 
M z 1.15 in the case of equation (37). 



Comparison of Results of Present Paper • 

v;ith Prandtl-Glauert Approximation 

The well-knovnn. Prandtl--Gla^v:vrt approxination Is 
based on the as siJinptioh of vani.::l'Lingl-^^ £=rn-f.ll disturb- 
ances to the main stream. Tho Prandtl^-Glauert velocity 
correction formula may be expressed as 



(**' — 

v.. 



^^1 



(57) 



where q - is vanisMngly small. The left-hand side 
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of this equation is actually the differential coeffi- 

clent -Ty — 7 — r-* evaluated at the main stream velocity 
d(q/^l)i 

q = q]_ .(or T = T^). An exact form of the Prandtl- 
Glauert approximation then is 



l(q/qi X 



T=T- 



(l - ) 



J72 



(58) 



The differential coefficient in equation (58) Is now 
evaluated for the various approximations treated in the 
present paper. 

For the arithmetic -mean approximation of the 
present paper given by equation (37) (y or B arbi- 
trary) , 



l(q/qi) 



i(q/qi)l 



T=T' 



2 



1 + (l - Mi2)(i + 
= 1 + I + I + ^ M^6 



128 p2 ^'^1 ^ 



(59) 
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For the Chaplygln or the Temple -Yarv^ood approxi- 
mation given by equation (41) (y = 1«4 or p = 2.5), 



"^(q/qi)e 



T=T 



^ - 4 ^1 
1 - 



Ml^ + |§ Mr + ••• (60) 



r J 



' For the von Karnan apiDroximatlon given by equa- 
tion (42) (y = -1 or p =*-^-), 



d 



d 



K-^l)c. 
(q/qi)l_ 



= (1 - T)l/2 



T=T- 



(l - Mi^) 



(61) 



For the geometric-mean approximation of the present 
paper given by equation (56) (y or p arbitrary). 



'^(q/qi); 




1 ■ 1 - T-^ . ■ 1 


fl/2 


d(q/qi)i 




1 - (2p + 1)t^ 






(62) 



Equation (62) Is indepencJent of the value of the 
adiabatic index y and includes the von Karman approxi- 
mationo Observe that the geometrlc-iaean approxiraation 
yields the Prandtl-Glauert result exactly, v\rhereas the 
arithmetic-mean approximation yields the Prandtl-Glauert 
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result Insofar as terms inclusive of Mj^ are concerned. 
The Ghaplygln or the Temple-Yarwood approximation con- 
tains the Prandtl-Glauert result only insofar as the M]_^ 
term is concerned. 



RESUME AND CONCLUDING RHiARKS 



1. Basic elementary solutions of the hodograph equa- 
tions have been employed to provide a basis for comparison, 
in the form of velocity correction formulas^ of corre- 
sponding compressible and incompressible flows. 

2. The velocity correction formulas obtained bj 
Chaplygin, by von Karmdin, and by Temple and Yarv/ood have 
been unified by means of these basic solutions and shown 
to be essentially equivalent • 

3. In the present paper tv/o types of approximations 
have been introduced by means of the basic elementary 
solutions, nam.ely, the ^*arithme tic-mean^* type and the 
"geometric-mean^^ type. These approximations include 
those obtained by Chaplygin, by von Karman, and by Temple 
and Yarwood. 

4* The approximations discussed in the present paper 
have been compared with the well-known results of Prandtl 
and Glauert. For this purpose, it has been emphasized 
that the Prandtl-Glauert result is valid for vanishingly 
small disturbances and, in a strict sense, is the slope 
term in a Taylor expansion in a qu.antlty which measures 
the disturbance. It v^as found that the arithmetic-mean 
type yields the Prandtl-Glauert result to a iiigher order 
of approximation than the Chaplygin or the Temple-Yarwood 
type and that the geometric-mean type contains the 
Prandtl-Glauert result exactly. The two types of approxi- 
m.ations introduced in the present paper then appear to be 
preferable to the others as a basis for extrapolation into 
the range of high stream Mach numbers and large disturb- 
ances to the m^ain stream. 

5. The results of the present paper have been ob- 
tained without consideration of any jjarticular boundary. 
The actual boundary problem of 'determ.ining the flow past 
a prescribed body is of a high order of difficulty and 
involves in general all the particular solutions of the 
hodograph equations. 
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6c The particular solutions discussed in the present 
paper are v/ell-behaved functions in both the subsonic and 
the supersonic regions. The hodograph equations give no 
reason, in general, to suppose that a. discontinuity neces- 
sarily occurs in the solution when local sound speed is 
attained. Rather,, it appears that the first breakdoTO of 
the 'solution is associated with the vanishing of the 
Jacobian of the transformation from^the physical to the 
hodograph variables. Indeed, von Karman has made an equi- 
valent suggestion in that the appearance of infinite accel- 
erations in the flow solution is a condition for flow dis- 
continuities. Interesting speculations on this matter are 
suggested by the results of the present paper since the 
"limiting^* curves discussed in the present paper are 
defined by a condition that is equivalent to the condition 
for infinite acceleration. The arithmetic-mean type of 
approximation thus yields a limiting value of the local 
Mach number M 1.15, and the geometric-Biean type of 
approximation yields a limiting value of the local Mach 
number M 1. The value M = 1 appears to be exact for 
vanishingly small disturbances! that is, local Mach 
nijmber M = stream Mach number ' = 1 (Prandtl-Glauert 

approximation) . However, for finite disturbances to the 

main flow due to the presence of a body in the fluid, 
infinite accelerations m.ay occur, for stream Mach numbers 
less than unity, in regions where the local Mach number 
is greater than unity. In this regard, the arithmetic- 
mean type of approxim.ation, considered as an extension of 
the Prandtl-Glauert relation to finite disturbances,- indi- 
cates the possibility of a mixed subsonic and supersonic 
flovv^ without discontinuities. It is important, however, 
to recognize that in general the limiting value of the 
local Mach number M is a function of shape parameters 
and is a result of the blending of many particular solu- 
tions of the hodograph flow equations according to the 
bouridary conditions. 
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x: 
o 


1.00 

.99750 
.99004 
.97769 

.93889 
.91271 
.89792 
; 882 13 
.86516 

.84700 
.83744 
.82756 
.81732 
.80664 
.79552 
.78386 
.77899 
.76880 
.75762 

1 


#^ 
to 
+ 


1.00 
.99760 
.99004 
.97769 
.96060 
.93890 
.91281 
.89817 
.88248 
.86578 . 

.84810 
.83889 
.82944 
.81975 
.80985 
.79969 
.78933 
.78511 
.77658 
.76791 

.75910 
.75015 
.74106 
! .73186 
.72261 
.71303 
.69869 
.68387 
.6739I& 
.62247 

.66865 
.51307 
.23869 
.05844 
.00496 
.00000 
.00000 

a 


• 


1.00 
.99749 
.98992 
.97711 
.95879 
.93463 
.90430 . 
.88671 
.86748 
.84656 

.82397 
.81202 
.79968 
.78691 
.77374 
.76016 
-.74617 
.74047 
.72888 
.71706 

.70498 
.69268 

.68017 
.66743 
.65449 
.64136 
.62130 
.60087 
.58706 
.61613 

.44356 
.37147 
.09096 
.00596 
-.00004 
.00000 
.00000 

0 


• 


1.00 
.99750 
•99016 
.97828 
.96241 
.94321 
.92140 
.90977 i 
.89775 
.88644 i 

.87295 
.86664 
.86032 
.85398 
.84764 
.84130- 
.83497 
.83245 
.82741 
.82237 

.81737 
.81239 
.80742 
.80260 
.79759 
.79273 
.78650 
.77834 
.77364 
.7507^ 

.72903 
.70866 
.62641 
.67164 
.63542 
.49392 
.47268 
.43162 




0 

-.00250 
-.01001 
-.02256 
-.04020 
-.06306 
-.09133 
-.10768 
-.12541 
-.14484 

-•16605 
-.17740 
-.18927 
-.20173 
-.21487 
-.22876 
-<i24353 
-.24976 
-.26292 
-.27757 


+ 


0 

-•00251 
-.01001 
-.02256 
-.04020 
-.06304 
-.09123 
-.10740 
-.12502 
-.14412 

-.16476 
-.17668 
-.18700 
-.19875 
-.21091 
-.22352 
-.23658 
-.24193 
-.25285 
-.26409 

-.27562 
-.28748 
-.29967 
-.31218 
-.32603 
-.33823 • 
' -.35869 
-.37999 
-.39464 
-.47406 

-.56448 
-.66734 
-1.43259 
-2 •83970 
-6.30024 
-16.98888 
-40.78337 
-OO 


to 


0 

-.00251 
-.01013 
-.02316 
-.04208 
-.06760 
-.10059 
-.12023 
-.14216 
-.16657 

-.19363 

-.20822 

-.22354 

-.23964 

-.25652 

-.27423 , 

-.29280 ' 

-.30047 

-.31624 

-.33261 

-.34958 
-.36718 
-.38542 
-.40432 
-.42390 
-.44417 
-.47594 
-.50938 
-.63263 
" -.66139 

-.81292 
-.99030 
-2.39742 
-6.12014 
-9.99177 
-31.27238 
-80.81740 
-OO 




0 

-•00250 
-.00989 
-.02196 
-.03831 
-•05847 
-.08186 
-.09457 
-.10787 
-.12167 

-.13588 
-.14313 
-.15045 
-.15786 
-.16530 
-.17281 
-.18036 
-•18339 
-.18946 
-.19556 

-.20166 
-•20778 
-.21391 
-.22003 
-.22616 
-.23228 
-.24144 
-.25059 
-•25665 
-.28673. 

-.31604 
-.34438 
-.46775 
-.65925 
-.62470 
-.70538 
-1.74934 
-.84019 




0 

.00200 
.00794 
.01768 
.03101 
.04762 
.06716 
.07792 
.08925 
.10112 

.11348 1 

.11982 j 

.12626 

.13279 

.13942 

.14612 

.15290 

.15663 

.16113 

.16667 

.17224 
.17786 
.18349 
1 * .18916 
.19485 
, .20056 
1 .20917 
.21782 
•22360 
.25262 

.28161 
. .31034 
.44444 
.55556 
.64286 
.76190 
.83333 
1.00 


• 
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TABLB 2.. VALnaa OP (q/qi)c, (q/tHu 



(q/«ii)e 



FOR r ■ 1*4 AND F0R-TARX0D8 VALDBS OP Mi 



NAIIONAL ADVISURY 
CUMMITUE FOR AERONAUTIC^ 





0.2 


0.3 


0.4 


0.5 


0.55 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


r 


0.00794 


0.01768 


0.03101 


0.04762 


0.05705 


0.06716 


0.08925 


0.X1M8 


0.19048 


0.16667 


0.19466 


0.88960 


Ml « 0,2 


(q/qi)c 


1.00 


1.49262 


1.97660 


2.44948 


2.68106 


8.90907 


3.36349 


3.78124 


4.19121 


4.58255 


4.96484 


6.30790 


(<i/qi)i 


Bq. (36) 


x.oo 


1.47471 


1.92120 


2.33332 


8.52518 


8.70706 


3.04051 


3.33365 


3.58793 


3!80599 


3.99122 


4.14788 


Eq. (36) 


1.00 


1.47331 


1.91444 


2.31266 


2.49242 


8.65746 


2.93870 


3.14735 


3.27693 


3.31937 


3.27590 


3*147.80 


Bq. (37) 


1.00 


1.47401 


1.91702 


2.32296 


2.60876 


8.66213 


2.98917 


3.23912 


3.42841 


3.56441 


3.61694 


3.61309 


(q/qi)i 


Eq. (37) 


1.00 


1.01263 


1.03066 


1.06446 


1.06668 


1.08461 


1.12168 


1.16797 


1.22249 


1.28926 


1.97026 


1.46907 


(Bq. fl8«)) 


0 


-1.17271 


-2.67803 


-4.39614 


-6.29383 


-6.19382 


-7.93514 


— 9.49190 


-10.76400 


-11.63383 


-12.07802 


-18.06448 




0 


•1.21250 


-2.82393 


•4.75600 


•6.81893 


-6.92214 


-9.89899 


-11.62814 


-14.00071 


-16.81967 


-18.61999 


-20.87000 


■l = 0.3 


(q/qi)o 




1.00 


1.32425 


1.64107 


1.79621 


1.94898 


8.24672 


2.69991 


2.80797 


9.07016 


3.31966 


3. 66618 


(q/qi)i 


Bq. (38) 


0.67810 


1.00 


1.30277 


1.66224 


1.71838 


1.88567 


* 8.06177 


2.26065 


8.48899 


2.58087 


2.70646 


2.81885 

— —f.. — ■ 


Bq. (36) 


0,67875 


1.00 


1.29942 


1.66971 


1,69172 


1.80374 


1.99463 


2.19626 


8.82355 


2.25904 


2.28358 


8.19666 


Bq. (37) 


0.67845 


1.00 


1.30110 


1. 67596 


1.70800 


1.81965 


2.02793 


2.19768 


2.92699 


2.41140 


2.46917 


2.46189 


(Q/<»l)e 




0.98753 


1.00 


1.01779 


1,04131 


1.05536 


1.07107 


1.10789 


1.15280 


1.20726 


1.27319 


1.96916 


1,46075 


(Bq. (18ft)) 


Bq. (37) 


0.53fl7S 


0 


-0.69286 


-1.48365 


-1.89680 


-2.81113 


-9.11250 


-3.62909 


-4.40996 


-4.81485 


-5*01604 


-5.00869 




0.55794 


0 


-0.74111 


-1.68984 


•8.11666 


-2 .62681 


-9.69898 


-4.7B981 


-5.68908 


-6.94746 


-7.86968 


-8.90671 




B 0.4 


(qAi)e 


0.50592 


0.75614 


1«00 


1.23924 


1.35640 


1.47174 


i;6to9 


1.9^900 


8.12040 


8.81840 


8.80874 


t.686a6 


(q/qi)i 


Bq. (85) 


0.52051 


0.76769 


1.00 


1.81461 


1.31438 


1.40903 


1.68860 


1,73519 


1.86754 


1.98106 


8.07744 


2.15866 


Bq. (»«) 


0.52235 


0.76967 


1.00 


1.80801 


1.30190 


1.38811 


1.69503 


1.64401 


1.71116 


1.73387 


1.71115 


1,64488 


Bq. (87) 


0.58U3 


0.76857 


1.00 


1.81124 


1.30813 


1.39853 


1.66868 


1.68897 


1.79r66 


1.85336 


1.68542 


1.8U94 


(«i/qi)c 
(q/qi)i 


Bq. (37) 


0.97025 


0.08253 


1.00 


1.02312 


1.03690 


1.06236 


1.06858 


1.18264 


1.18614 


1.26092 


1.32964 


1.48940 


(Bq. il8ft)) 


0.72811 


0.40930 


0 


-0.46710 


-0.71120 


.0.95589 


-1.48980 


-1.86262 


-2.19669 


-2.43494 


-2.56481 


-8.64929 


(bb. (lab)) 


0.76643 


0.43714 


0 


-0.62420 


-0.81188 


•1.11259 


-1.74168 


•9.88866 


-9.09429 


-9.66205 


•4.25988 


•4.81760 






Ml e 0.8 






(q/qi)o 


0.40825 


0.60936 


0.60695 


1.00 


1.09454 


1.18763 


1.96906 


1.64970 


1.71106 


1.87084 


2.02281 


2.16694 


(q/qi)i 


Bq. (35) 


0.42867 


0.66202 


0.82336 


1.00 


1.08222 


1.16017 


1.90907 


1.42871 


1.59769 


1.63115 


1.71061 


1.77737 


Bq. (36) 


0.43240 


0.63706 


0.62781 


1.00 


1,07773 


1.14909 


1.27069 


1.36092 


1.41662 


1.43531 


1.41651 


1.96110 


Bq. (37) 


0.43049 


0.63454 


0.82560 


1.00 


1.07998 


1.15462 


1.88679 


1.39440 


1.47588 


1.53012 


1.68668 


1.55699 


(q/qi)c 
TqTqJTJ 


Bq, (37) 


0.94634 


0.96032 


0.97741 


1,00 


1.01348 


1.02859 


1.06999 


, 1.10707 


1.15935 


1.22867 


1.89941 


1.90918 


(Bq. (18*)) 


0.81468 


0.59736 


0.31838 


0 


-0.16636 


-0.33316 


-0.66669 


.0.94435 


-1.17822 


-1*. 94127 


-1.42907 


-1.41924 


(BQ.^tSb)) 


0.BT771 


0,65366 


0.36646 


0 


-0.19554 


.0.40000 


-0.68766 


.1.26764 


-1,70667 


-2.13349 


-2.53960 


•2.91900 






Hi « 0.55 








(q/qilo 


0.37299 


0.6567S 


0.73725 


.0.91363 


1.00 


1.08604 


1.25081 


1.41036 


1.66326 


1.70929 


1.84810 


1.97978 


(q/qi)i 


Bq. (36) 


0.39601 


0.58399 


0.76083 


0.92404 


1.00 


.1.07202' 


1.20408 


1.32017 


1.42086 


1.60722 


1.68057 


1,64235 


Kq. (36) 


0.40122 


0.69110 


0.76811 


0.98788 


1.00 


1.06621 


1.17906 


1.26278 


1.31434 


1.93176 


1 .91495 


1.26294 


Bq. (37) 


0.39861 


0.68753 


0.76446 


0.98595 


1.00 


1.06911 


1.19160 


1.29114 


1.36657 


1.41680 


1.44193 


1 .44019 . 


(q/qi)c 
Tq7qp7 


Bq. (37) 


0.93573 


0.94786 


0.96441 


0.98669 


i.od 


1.01490 


1.0497JB 


1.09234 


1.14393 


1.20640 


1.26222 


1.37467 


Cp,o 
(Bq. (18«)) 


0.84111 


0.65481 


0.41560 


0.14262 


0 


-0.14300 


-0.41967 


-0,66704 


-0.66781 


-1.00792 


-1.07743 


-1.07418 


(Bq. (IBb)) 


0.91635 


0.72685 


0.47249 


0.16T67 


0 


-0,17497 


•0.54078 


-0.91791 


-1.29899 


-1.66828 


-8.00690 


•8.99068 
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NAIIUNAL AOVISUKV 

TABLE 8.- Contlnu.4 COMMIIUt FOU AtRONAUllCS 



H 


0.4 


0.6 


0.65 


0,7 


0,76 


0,8 


0.85 


0,9 


1.0 


1.1 


. 1.2 


1.3 


T 


0,03101 


0.06716 


0,07798 


0.08925 


0.10112 


0.11348 


0.12626 


0.13942 


0.16667 


0,19485 


0.88360 


0.86868 


= 0,6 




0.67947 


1.00 


1.07707 


1.15277 


1.22704 


1,29981 


1.37107 


1.44093 


1.57527 


1.70324 


1,82460 


1.93939 




Kq. (36) 


0.70971 


1.00 


1.06348 


1.12518 


1.17914 


1.23146 


1.28018 


1.32659 


1,40596 


1.47436 


1.53201 


1.68014 


(q/qi)i 


Bq. (S6) 


0.72041 


1.00 


1.06612 


1.10583 


1.14869 


1,18435 


1.21245 


1.23289 


1,24908 


1.23272 


1.18461 


1.10691 


Bq. (37) 


0.71504 


1.00 


1.05979 


1.11447 


1.16562 


1,20767 


1.24585 


1.27641 


1.32521 


1.34615 


1.34708 


1.32253 


(q/qi)c 
(q/qiii 


Bq. (37) 


0.95025 


1.00 


1,01631 


1.03437 


1.05432 


1.07630 


1.10051 


1.12713 


1,18869 


1.26339 


1.55449 


1.46642 


(Bq. (18*)) 


0,48872 


0 


-0,12315 


-0 .24204 


.0.35448 


•0.45847 


-0.55214 


-0.63433 


-0.75618 


-0.81751 


-0.81462 


-0,74909 


<=p,Ki 


0.6649S 


0 


-0.16786 


-0.31929 


«iO«48910 


.0.64774 


.0.61856 


.0.97699 


.1.89437 


-1.59768 


-1.88099 


*»8. 14151 
















« 0,65 












(q/qi)c 


0.650B4 


0.92843 


1,00 


1.07039 


1.13924 


1,20681 




1.33765 


1.46255 


1.56158 


1.69405 


1,60060 


(q/qi)i 


Bq. (55) 


0.66734 


0.94030 


1.00 


1.06615 


1.10878 


1.15797 


1.20577 


1.24630 


1.32204 


1,38638 


1.44057 


1.48582 


Bq. (36) 


0.68212 


0.94685 


1.00 


1.04T08 


1.08766 


1,12142 


1.14802 


1,16723 


1.18271 


1,16725 


1,12158 


1.04808 


Bq. (37) 


0.67469 


0.94356 


1.00 


1.05169 


1.09816 


1.13955 


1.17555 


1.20612 


1.25044 


1,27209 


1,27110 


1.24789 


(q/qi)e 
(q/qi)i 


Bq. (57) 


0.95501 


0,98396 


1.00 


1.01778 


1.03741 


1.05904 


1.08286 


1.10906 


1,16965 


1.24314 


1.33274 


1.44292 


(Bq, (IBa)) 


0.54479 


0.10969 


0 


-0.10684 


-0.20596 


-0.29853 


-0.38192 


.0.46473 


-0.56360 


-0.61821 


-0.61570 


-0.65723 


(BQ. (18b)) 


0.64074 


0.14O0e 


0 


.0.14863' 


•0.28e62 


-0.43476 


.0.68086 


•0.78671 


-1.00882 


.1.27767 


.1.68915 


-1.76084 










«3 


= 0,7 






(q/qi)e 


0.58042 


0.86747 


0.93453 


1.00 


1,06442 


1.12766 


1.18936 


1.24980 


1.56680 


1.47762 


1.68280 


1.68235 


(q/qi)i 


Bq. (35) 


0.65186 


0.69032 


0.94685 


1.00 


1.04965 


1.09641 


1.13978 


1.18008 


1.25175 


1.51266 


1.36598 


1.40681 


Bq. (35) 


0.65146 


0.90429 


0.95504 


1.00 


1.03875 


1.07099 


1.09640 


1.11475 


1.12955 


1.11476 


1,07114 


a. 00096 


Bq. (57) 


0.64159 


0.69728 


0.96094 


1,00 


1.04429 


1.08362 


1.11788 


1.14695 


1.16909 


1.2()969 


1.20874 


1.18666 


(q/qi)c 
(q/^iii 


Bq. (57) 


0.91869 


0.96678 


0,98255 


1,00 


1.01926 


1,04055 


1,06394 


1.08967 


1.14920 


1.22140 


1,30946 


1,41772 


(Bq. (le*)) 


0.58856 


0.19489 


0.09571 


0 


.0.09054 


-0.17423 


-0.24966 


.0.31549 


-0.41394 


-0,46555 


-0,46105 


-0.40816 


Cp.Ml 
(Bq. (leb)) 


0.70580 


0.25516 


0.12904 


0 


.0.15082 


-0.26254 


-0.39997 


-0.52440 


-0.77907 


-1,02140 


-1,24761 


-1.46569 




Ml = 0,76 






■ (q/qi)c 


0.56374 


0.81497 


0.87778 


0.93947 


1.00 


1,05931 


1.11757 


1.17416 


1.28380 


1,38810 


1,46701 


1,58054 


(q/qi)i 


Bq. (35) 


0.60187 


0.84807 


0.90191 


0.95253 


1,00 


1,04436 


1.08566 


1.18404 


1.19255 


1.25056 


1,29926 


1.34006 


Bq. (56) 


0.62715 


0.87056 


0.91942 


0,96868 


1.00 


1,03105 


1.05560 


1.07317 


1,08740 


1.07317 


1.03120 


0.96562 


Bq. (37) 


0.61459 


0.85924 


0.91062 


0.96759 


1,00 


1.03766 


1.07047 


1.09831 


1,15668 


1.15840 


1.15749 


1.13636 


(q/qi)c 
(q/qi)i 


Eq. (37) 


0.90126 


0,94848 


0.96394 


0.98108 


1.00 


1.02084 


1,04381 


1,06906 


1,12745 


1.19629 


1.28468 


1,39088 


(Eq. (lea)) 


0.62252 


0.26171 


0.17077 


0,08302 


0 


-0.07678 


-0,14591 


.0.20628- 


-0.29659 


-0.34189 


-0.33976 


-0.29151 


Cp,Ml 
(BQ. (18b)) 


0.763%! 


0.35197 


0.23700 


0.11957 


0 


-0.12006 


-0,23997 


-0.35895 


-0.59124 


-0.61227 


-1.01672 


-1.20653 


Ml =f 0.8 










(q/qi)c 


0.52274 


0.76934 


0.82864 


0,68687 


0.94402 


1.00 


1,05482 


1.10842 


1.21195 


1.31036 


1.40376 


1.49204 


(Q/qi)i 


Bq. (35) 


0.57631 


0.61205 


0.86359 


0.91207 


0.95755 


1.00 


1.03956 


1.07629 


1.14171 


1.19726 


1.24406 


'1.28512 


Bq. (36) 


0.60827 


0.64435 


0.89172 


0.93571 


0.96990 


1.00 


1.02372 


1.04085 


1.05467 


1.04086 


1,00014 


0.95460 


b4. (37) 


0.59209 


0.82606 


0.87756 


0.92265 


0.96570 


1.00 


1.03162 


1,05845 


1.09734 


1.11635 


1.11545 


1.09610 


(q/qi)e 

(q/qi)i 




0,88287 


0.92910 


0.94426 


0.96105 


0.97968 


1.00 


1.02249 


1.04723 


1.10443 


1.17383 


1^25646 


1.56247 


Cp,o 
(Bq. (18a)) 


iSq. (57) 


0.64943 


0.31433 


0.22989 


0.14858 


0.07126 


0 


-0.06424 


-0.12027 


-0.20416 


-0.24619 


-0.24423 


-0.19984 


<^P.Ml 
(Bq. (IBb)) 




0.81520 


0.43547 


0.32940 


0.22085 


0.11074 


0 


-0,11065 


•0.22036 


-0.43464 


-0.63857 


-0.62904 


-1.00411 



o 
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TABLE Concluded 



NAIIUNAL ADVISURY 
CUMMlTIEe FOR AERONAUTICS 





0.4 


0.6 


0.825 


0.85 


0.675 


0.9 


0.925 


0.96 


1.0 


1.1 


1.2 


1.3 


r 


0.03101 


0.06716 


0.11982 


0.12626 


0.13279 


0.13942 


0.14612 


0.15563 


0.16667 


0.19485 


0.22360 


0.26262 










Ml 


=0.825 














(q/qi)e 


0.50872 


0.74871 


1.00 


1.02663 


1.05276 


1.07870 


1.10434 


1.13972 


1.17942 


1.27624 


1.36610 


1.45203 




Bq. (56) 


0.56494 


0,79602 


1.00 


1.01905 


1.03738 


1.06506 


1.07205 


1.09476 


1.11917 


1.17563 


1.21950 


1.25781 




Bq.. (36) 


0.60068 


0 .83360 


1.00 


1.01093 


1,02021 


1.02785 


1.03380 


1.03929 


1.04147 


1,02784 


0.98764 


0.92292 




Eq. (57) 


0.58253 


0.81469 


1.00 


1.01497 


1.02674 


1.04136 


1.05273 


1.06665 


1.07962 


1,09833 


1.09746 


1.07744 


^ {q/qi)c 




0.67329 


0.91901 


1.00 


1.01139 


1.02335 


1 * 03586 


1.04902 


1.06850 


1 .09244 


1 .16107 


1 .24476 


1 .34767 


Cp,o 
(BQ. ilM) 


Bq. (57) 


0.66066 


0.33628 


0 


-0.03016 


-0.05831 


-0.06443 


-0.10824 


-0.13774 


-0.16558 


-0.20633 


-0.20442 


-0.16088 


(&. (18b)) 




0.859S0 


0,47330 


0 


-O, 05329 


-0.10631 


-0.15910 


-0.21149 


-0.28409 


-0.36571 


-0.56234 


-0.74599 


-0.91467 








Ml 


= 0.85 












(q/qi)c 


0»49S&8 


0.72936 


0.97416 


1.00 


1.02555 


1.05062 


1.0^580 


1.11026 


1.14894 


1.24228 


1.33080 


1.41451 


(q/qi)i 


Bq. (55) 


0.55458 


0.78114 


0.96132 


1.00 


1.01799 


1.03532 


1.05202 


1.07430 


1.09826 


1.15171 


1.19672 


1.23432 


Eq, (56) 


0.59418 


0.82476 


0.98920 


' 1.00 


1,00917 


1.01673 


1,02264 


1.02805 


1.03023 


1.01673 


0.97697 


0.«1295 


Bq. (57) 


0.67594 


0.80267 


0.98526 


1.00 


1 ..01357 


1,02599 


1,03721 


1.05091 


1.06371 


1.08212 


1.08128 


1.06156 


(q/qi)c 
(q/qi)i 


Eq. (37) 


0.86347 


0.90867 


0.98873 


1.00 


1*01182 


1.02420 


1.03721 


1.05647 


1.08013 


1,14801 


1.23076 


1.33249 


Cp,o 
(Bq. {18a)) 


0.67069 


0.35572 


0.02926 


0 


-0,02732 


-0.05266 


-0.07580 


-0.10441 


-0.13148 


-0.17098 


-0.16917 


-0.12689 


Cp,Mi 
(Bq. (18b)) 




0.86288 


0.50895 


0.05147 


0 


«0. 05129 


-0.10226 


-0.15294 


-0.22305 


-0.30199 


-0.49202 


-0.66952 


-0.83274 












= 0.876 










(q/qi)c 


0.48322 


0.71118 


0.94988 


0.97508 


1.00 


1.02463 


1.04898 


1.08259 


1.12031 


1.21133 


1,29763 


1.37926 


(q/qi)i 


Bq. (55) 


0.54458 


0.76734 


0.96396 


0.98232 


1.00 


1.01703 


1.03340 


1.06530 


1.07886 


1.13134 


1.17555 


1.21248 


Bq. (56) 


0.68877 


0.81727 


0.98020 


0.99091 


1,00 


1.00748 


1.01333 


1.01870 


1.02085 


1.00749 


•0.96807 


0.90464 


Eq. (57) 


0.56626 


0.79192 


0.97207 


0.98661 


1.00. 


1.01225 


1.02331 


1.05684 


1.04947 


1.06764 


1.06678 


1.04733 


(q/qi)c 
(q/qi)i 


Bq. (57) 


0.85535 


0.89805 


0.97717 


0.98831 


1.00 


1.01223 


1.02509 


1.04412 


1.06750 


1.13459 


1.21640 


1.31693 


Cp,o 
(Bq. (18a)) 


0.67935 


0.37286 


0,05608 


0.02660 


0 


-0.02465 


-0.04716 


-0.07504 


-0.10139 


-0.13986 


-0.13602 


-0.09690 


Cp.Mi 
(Ba. (18b)) 


0.88564 


0.54277 


0.09960 


0.04971 


0 


-0.04943 


-0.09846 


-0.16642 


-0.24286 


-0.42701 


-0.69898 


-0.75712 






Ml = 0.9 








(q/qi)c 


0.47160 


0.69409 


0.92705 


0.95164 


0.97595 


1.00 


1.02376 


1.05657 


1.09338 


1.18220 


1.26644 


1.34609. 


<q/qi)i 


Eq. (55) 


0.55546 


0.75449 


0.94783 


0.96587 


0.98325 


1.00 


1.01610 


1.03764 


1.06079 


1.11240 


1.15B88 


1.19218 


Eq. (36) 


0.58459 


0.81121 


0.97291 


0.98555 


0.99257 


1.00 


1.00579 


1.01113 


1.01327 


. 1.00000 


0.96069 


0.89792 


Eq, (37) 


0.56939 


0.78233 


0.96029 


0.97466 


0.98788 


. 1.00 


1.01091 


1,02429 


1.03676 


1.05470 


1.05387 


1.03463 


(q/qi)c 
(q/qi)i 


Eq.. (37) 


0.84306 


0.88721 


0.96539 


0.97638 


0.98792 


1.00 


1.01271 


1.03151 


1.05461 


1.12089 


1.20170 


1.30104 


Cp,o 
(Bq. (18a)) 


0.68708 


0.38796 


0.07784 


0.05004 


0.02409 


0 


-0.02194 


-0.04917 


-0.07487 


-0.11239 


-0.11064 


-0.07046 


Cp.Mi 
(Bq. (18b)) 




0.90787 


0.67490 


0.14460 


0.09617 


0.04799 


0 


-0.04762 


-0.11365 


-0.18788 


-0,36660 


-0.53362 


-0.68711 






Ml = 0.925 






<q/qi)c 


0.46066 


0.67798 


0.90552 


0.92955 


0.95331 


0.97678 


1.00 


1,03205 


1.06800 


1.15476 


1.23704 


1.31486 


(q/qi)i 


Bq. (55) 


0.52697 


0.74253 


0.93279 


0.95057 


0.96768 


0.90415 


1.00 


1.02119 


1.04397 


1.09476 


1.13756 


1.17329 


Eq. (36) 


0.58103 


0.80654 


0.96729 


0.97787 


0.98686 


0.99423 


1.00 


1,00532 


1.00743 


0.99424 


0.95534 


0.89276 


Eq. (37) 


0.56335 


0.77387 


0.94991 


0.96413 


0.97721 


0.98919 


1.00 


1.01323 


1.02556 


1.04332 


1.04249 


1.02347 


(q/qi)c 
(q/qi)! 




0.83249 


0.87609 


0.95327 


0.96413 


0.97554 


0.98745 


1.00 


1.01857 


1.04138 


1.10681 


1.18662 


1.28471 


Cp,o 
(Eq. (18a)) 


Bq. (37) 


0.69380 


0.40113 


0.09767 


0,07045 


0;04506 


0.02150 


0 


-0.02664 


-0.05177 


-0.08852 


-0.08679 


-0.04749 


Cp.Mi 
(Eq. (18b)) 




0.92963 


0.60574 


0.18708 


0.13991 


0.09301 


0.04633 


0 


-0.06425 


-0.13642 


-0.31038 


-0.47285 


-0.62222 



••. :•• : • •*/ : :V :*: 
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mBLE 3.- VALUES OP P(t) FOR SEVERAL VALUES OP y 



T) = 1 - (gP + 1)T _ 1 - 





F 




Y = 1.4 


Y = 1 


Y = 2 


Y z: 00 


M 


(P = 2.5) 




(P - 1) 


(p->o) 




Adiabatic 


Isothermal 


Hydraulic 


Limiting 








analogy 


incompressible 






(1) 


(2) 


0 


1 .00 


1 00 


1 on 


1 DO 


.2 








» youu 


.4 


.98328 


.98575 


.97977 


,8400 


.6 


.90606 


.91733 


.89113 


.6400 


.65 


.86634 


.88113 


.84726 


. .5775 


.70 


. 81394 


. 83248 


.79050 


.5100 


.75 


.74558 


.76783 


.71822 


.4375 


.80 


.65738 


. 68273 


.62728 


.3600 


.85 


.54489 


.57153 


.51421 


.2775 


.90 


.40258 


.42710 


-.37504 


. 1900 


.95 


.22355 


.24041 


.20534 


.0975 


1.00 


0 


0 


0 


0 


1.05 


-.27752 


-.30870 


-.24667 


-.1025 


1.10 


-.62069 


-.70423 


-.54102 


-.2100 


1.20 


-1.55960 


-1.85711 


-1.30158 


-.4400 


1.30 


-2.95915 


-3.73944 


-2.34862 


-.6900 


1.50 


-8.01227 


-11.8597 


-5.64466 


-1.2500 


2.00- 


-56.6884 


-163.79 


-26.9981 


-3 . 0000 



1 Y = 1, p =r (1 - M^)e^^ 



2y = ", F=1-m' 
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Fig. lb 
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Figs. 2a. b 
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Fig. 3 
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Fig, 6 




